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Response functions for spin-density-wave (SDW) and d-wave singlet superconductivity (dSC) 
in quasi-one-dimensional (QID) electron systems are calculated by a renormalization group 
technique. It is shown that the response functions for both SDW, Xs, and dSC, Xd, are enhanced 
by interchain hopping, t^, i.e., by quasi-one dimensionality. When the Fermi surface deviates 
from perfect nesting, Xs saturates below the energy scale of imperfectness of the nesting, while 
Xd is hardly affected. Consequently, the superconducting correlation increases even away from 
the SDW phase. This gives a possible interpretation of the recent experimental results of QID 
organic conductor (TMTTF)2SbF6, where Tc increases even away from the SDW phase. 

KEYWORDS: superconductivity, quasi-one-dimension, renormalization group, dimensionality, spin fluctu- 
ation, superconducting fluctuation, organic superconductor, (TMTTF)2SbF6 



When a superconducting phase appears next to some 
ordered phase, it is reasonable to consider that the mech- 
anism of the superconductivity is due to the fluctuation 
of the ordered phase. Since the amplitude of fluctua- 
tion is reduced as the system is away from the ordered 
phase, we expect that the transition temperature of su- 
perconductivity, Tc, is also reduced monotonically. For 
example, in quasi-one-dimensional (QID) organic con- 
ductors TMTSF- and TMTTF-salts, its superconduct- 
ing phase appears next to the spin-density-wave (SDW) 
phase, and Tc decreases monotonically as the system 
leaves away from the SDW phase by applying pressure.^ 
This situation have been considered to be a typical ex- 
ample of the spin-fluctuation mediated superconductiv- 
ity.^"^ The pressure-temperature (P-T) phase diagram 
in these salts have been interpreted by the simple spin- 
fluctuation mechanism.'^ 

However, a recent careful measurement of 
(TMTTF)2SbF6 shows that Tc increases by pres- 
sure even away from the SDW phase. ^ This behavior 
of Tc cannot be understood from the above simple 
spin-fluctuation mechanism. It suggests a possibility 
that some other unknown effects enhance Tc. 

Our idea here is to take account of the dimensionality 
change due to the pressure as a key to understand this 
mysterious behavior. So far, the effect of pressure has 
been considered to be just to destroy nesting properties 
of the Fermi surface and to reduce the spin-fluctuation. 
However, the pressure simultaneously increases the di- 
mensionality through the increase of interchain hopping, 
t± . The change of the dimensionality possibly affects the 
superconducting correlation. A pioneering work of this 
subject was done by Suzumura and Fukuyama,^ who 
showed Tc of full gapped (s- or p-wave) pairing is an 
increasing function of t^. Though their argument is ap- 
plicable to the system with attractive interaction, it can- 
not be directly applied to the repulsive cases. There are 
some issues to be taken into account: interchain inter- 
action, anisotropic pairing, such as d- or /-wave, and 



competition between superconductivity and SDW. These 
points have been difficult to consider theoretically so far. 
However, due to the improvement of the renormalization 
group (RG) technique, we are now able to study these 
problems in a reliable way. 

In this letter, we investigate spin- and 
superconducting-fluctuation in QID electron sys- 
tems by changing the dimensionality and the nesting 
conditions with the newly developed RG technique for 
QID systems called A^-chain RG.* In particular, we 
clarify the following points: 1) is it possible for Tc to keep 
increasing even away from the SDW phase?, 2) where 
is the maximum of Tc as a function of pressure or of 
dimensionality?, 3) what are the roles of dimensionality 
and nesting property? In order to discuss the effect of 
dimensionality, it is of the prime importance to deal 
with ID fluctuation adequately, which results from the 
quantum interference between Cooper (particle-particle) 
and Peierls (particle-hole) channels. Since this is difficult 
in RPA or FLEX type approximation,'^''' we employ the 
RG method, which has been approved for taking into 
account the ID fluctuation.^'^" The obtained results 
indicate that the superconducting correlation essentially 
increases for small tj_, and can keep increasing even 
when the spin fluctuation is suppressed. 

We study the Hubbard model in weakly-coupled A^ 
chains with length L. In particular, we investigate this 
model at quarter-filling, which models TMTSF- and 
TMTTF-salts. The present results, however, is valid 
for other fillings except for the half-filling. For the 
convenience of the RG formulation, we express the 
Hamiltonian in the form H — Hq + Hi with Hq = 
Ep,k,<T^p(k)c];,k,^Cp,k,<T and 
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+ 94pPp{<l)Pp{-<i) + giaSpiqjS^i-q)]. (1) 

Here the coupling constants are given by gp — g^p = 
—g^r = —gAa = U for the Hubbard model. Hi = 
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U^^niiTiii, and Cp^k.o- (Cpku) annihilation (cre- 

ation) operator with spin a close to the Fermi surface 
of the right {p = 1) and left {p — —1) branches. The 
charge and spin density operators of the branch p are 
given by pp(q) = (1/2) Ek,^ i^.^^ ,^Cp,k,a and Sp(q) = 

(1/2) I]k,a./3 4-k+q,a'^"'''^"-'^-'^- ^^^^ Hamiltonian 



eq. (1), 8fc_F-Umklapp process is discarded. In the low 
pressure region of TMTTF-salts, this Umklapp process 
is relevant and leads to the charge order transition. 
In the high pressure region, on the other hand, the ef- 
fect of the 8fci?-Umklapp process is naively expected to 
be small. Thus we investigate the above Hamiltonian for 
the first step. 

We use the quasiparticle dispersion given by 

^p(k) = vp{pk — kp) — 2t± cos k± — 2t±2 cos 2k±, (2) 

where vp{— \/2t) is the Fermi velocity at quarter-filling 
with t being the intrachain hopping, kp being the Fermi 
wave number in the limit of t±^ — 0. and t\ji are 
the nearest and the next-nearest interchain hopping, re- 
spectively. When > and t\_2 = 0, the Fermi sur- 
face is perfectly nested with the nesting vector Qo = 
(±2^;^, ±7r) even though the Fermi surface is warped. On 
the other hand, when t±,t±2 > 0, the nesting becomes 
imperfect. These two parameters represent different ef- 
fects. t± corresponds to the change of dimensionality, 
where the crossover from ID to QID occurs by keeping 
the nesting properties, while t±2 breaks the nesting con- 
dition. In the real compounds, we can change t±_ and 
t±2 by applying pressure P. The relation between them 
depends on the various conditions of the experiments. 

We apply the A^-chain RG to the above Hamiltonian. 
This technique has a strong point that it can deal with 
the ID fluctuation effect even in the region T <t±, where 
the system is not a non- Fermi liquid state within the one- 
loop approximation. The details of this RG method are 
described in ref. 8. The response functions are calculated 
from 
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where Iqp = tanh[_E/4r], is pair vertex part, and 
^ s for SDW and = d for dSC. The initial band- 
width, Eq = 2vpkp, is renormalized as i? = Eoe~^. The 
pair vertex part of the SDW is renormalized as 
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Here, (^^(fci^, A:2±, fcax; ^) is a renormalized four-point 
vertex at step £, where fcix, —k2± and fcax are perpen- 
dicular components of incoming momenta ki,— k2, and 
outgoing momenta ka , respectively. Note that k4 is given 
by — k4 = ka + k2 — ki. Generally, marginal interactions 
in the low-energy region depend only on polar coordi- 
nates (or in this case, ki±) of the wave vectors. This is 
justified by the zero-order scaling analysis. Thus the 
dependence of on momenta parallel to the chain is 
irrelevant. The initial value for the vertex is given by 
G^,{£ = 0) = g^/iTVp, where gs = (-.gp + go-)/4 for SDW. 
For superconductivity, we consider partial components 




le-05 lc-04 0.001 0.01 0.1 

Temperature T / 1 

Fig. 1. Response functions of SDW (solid lines) and dSC (dashed 
lines) as a function of temperature for the perfectly nested Fermi 
surface tx2 = with U = 3t,N = 15, = 0,0.01,0.1. Inset 

shows the response functions for the imperfect nesting case with 
tx2 = 0-U±. 
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Fig. 2. Energy flows of the coupling Go-(fc±i, ^xa) for (a) 
tj_ = O.Olt and (b) tj^ = O.lt with 4^2 = 0, U = 3.0t at T = 
lD-6t. 



of the pair vertex. The pair vertex part of dSC is renor 
malized as 

d 



d£ 



\iizdi£)^~adlcpii), 



(5) 



where fld = [A/N'^) J2k^,k'^ cosk±Gsin{k±,k±,k'j_) cosk'j_ 
and gsin — {gp — iga)/^- (For the details of pairing 
symmetry, see ref. 19.) 

Figure 1 shows the temperature dependence of the re- 
sponse functions, XmC^)' f^'^ = 0.0,0.01,0.1 with 
U — 3t,N — 15. First we concentrate on the case 
where the Fermi surface is perfectly- nested, i.e., t±2 = 0. 
For purely ID systems {t± = 0), it is well known that 
Xs{T) shows a power- law behavior as Xs{T) oc T^"" with 
ao ~ U/2iTVp. In the case of QID systems (<_l > 0), 
Xs{T) is enhanced from that of ID, keeping apparently 
the power-law behavior as shown in Fig. 1. The power a 
of QID Xs{T) increases as t± increases and saturates for 
relatively large t± as shown in the inset of Fig. 3. The 
response function for dSC, Xd{T), is also enhanced by t± 
and tends to diverge at low temperatures, though it does 
not show any instability at = due to = 0.^^ 
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Fig. 3. Response function of QID normalized by that of ID, 
X^/Xmi ^ ^ function of t± for tj_2 = 0, U = 3t and T = O.Olt. 
The inset shows the power of Xs(T) (ex T"") for U/t = 2, 3. 



The instability of dSC originates from the k±- 
dependence of the coupihngs. Figures 2 show some exam- 
ples of energy flow of —Ga{ki±, ki±,k3±) for several sets 
of (fci^, k^±), which corresponds to the Cooper pair scat- 
tering (ki, — ki) (ks, — ka), with t±2 = and U = 3t. 
In the high energy region E > E^^ ^ 0{t±), every flow is 
identical with that of ID. In this region, interference be- 
tween Cooper (C) and Peierls (P) channels is dominant 
even though > 0, so that the couplings flow as those 
of ID. For E < Ex, on the other hand, the CP interfer- 
ence becomes imperfect, namely, C and P do not cancel 
each other in some types of coupHngs.^ As a result, each 
couplings flows differently and G^{ki±, k2±, fcsi) has k±- 
dependence. The amplitude of the fc_L-dependence of the 
couplings becomes larger as t± becomes larger. This k±- 
dependence causes the increase of Xd- In other words, 
the increase of Xd is due to the breakdown of the perfect 
CP interference. This is a characteristic behavior of the 
QID systems. In two or three dimensions, there are no 
perfect CP intereference in high-energy region, and thus 
the increase of Xd in the present context does not occur. 
Actually, for the isotropic square-lattice, the RG method 
yields the monotonically decrease of T^.^^'^"^'^^ 

Both Xs and Xd increase for small t± as we have seen, 
but there are difference for large t± . We plot the response 
functions normalized by that of ID, x^, as a function of 
t± for the perfect nesting case {t±2 = 0) at T — O.Olt in 
Fig. 3. For small t± (< 0.25t), both Xfi/Xfj. are propor- 
tional to t\. For large t± (> 0.25t), on the other hand, 
X/j/Xa* behaves differently: Xd/Xd saturates and sHghtly 
decreases, whereas Xs/Xs keeps increasing. These differ- 
ent properties of Xi^lXti against ij_ can be understood 
intuitively as follows. For SDW, the exchange coupling 
between chains stabilizes the long-range SDW order, 
and enhances x&- This process is a single-particle process 
as is shown in Fig. 4 (a), so that Jg ^ i\/U . Here, is 
defined as = (x^/Xm ~ 1) /Xa" which can be regarded 
as an effective interchain-coupling. For dSC, on the other 
hand, the interchain Josephson coupling stabilizes the 
dSC order and enhances Xd- This process is a pair pro- 
cess as is shown in Fig. 4 (b). In this case, Jd ^ t^/A^, 
where is the energy to break the local intrachain d- 
wave pair. When t±^ becomes larger {t± > 0.25), elec- 




Fig. 4. Illustration of the interchain process for (a) SDW and 
(b) dSC. The long-range SDW order is stabilized by the inter- 
chain exchange coupling, which is the single-particle process. The 
long-range dSC order is stabilized by the interchain Josephson 
coupling, which is the pair process. The coherent pair hopping 
is suppressed for large t± . 








SDW 




_^40 — 






dSC 




^7'"^ 




--120--5!; ..^ ^ 















Interchain hopping 1^/ 1 



Interchain hopping 1^/ 1 



(b) 




0.1 0.2 0.3 0.4 05 

Interchain hopping t^/ 1 



Fig. 5. Contour plot of TTVpXii{T,t±) with different nesting con- 
ditions: (a) tj_2 = 0, (b) tj_2 = — cos fcf'/4f sin^ /c^, (c) 
tx2/t (X [t^/t)'^. The insets of (b) is the linear plot of the low- 
temperature region. 



trons hop between chains not as a pair but individually, 
so that the coherent pair hopping is suppressed. This is 
why only Xdlxd saturates for large whereas Xs/Xs 
keeps increasing due to the single-particle process. From 
this viewpoint, both intra- and inter-chain Cooper pair 
forms impartially for t±_ > 0.25. Since both the intrachain 
(ID) spin-fluctuation and the interchain Josephson cou- 
pling help to form pairs, Xd have a maximum at around 
t±_ - 0.25. 

Next, we show the contour plot of TrvFXiJ.iT,t±) in 
Fig. 5. We plot three different situation of the nesting 
property, which correspond to different situation of the 
pressure effect. First, let us look into the role of the di- 
mensionality exclusively by assuming t±2 = in Fig. 5 
(a). In this case, the Fermi surface is perfectly nested, 
and thus Xs is always larger than Xd and keeps increas- 
ing by t± . Note that the contour lines of Xd has a broad 
maximum and turns to decrease at around t± ~ 0.25t 
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even though the Fermi surface is perfectly nested and 
the spin-fluctuation is not weakened. This decrease is 
due to the individual interchain hopping of the electrons 
mentioned above/ 

In the actual compounds, the pressure also has a role 
of increasing t±2 and destroying the nesting property. 
For example, when we linearlize the tight-binding form 
for the anisotropic square lattice, ^(k) = —2tcosk — 
2t±cosk±^ — II, we obtain t±2 = — 1\ cos kp/ At sin^ kp- 
Wc plot x,AT-t±) with this condition in Fig. 5 (b). The 
contour lines of Xs suddenly fall off when they cross the 
line T = t±2, since Xs{T) is cut off at T = t±2 as is shown 
in the inset of Fig. 1. On the other hand, Xd not af- 
fected so much by the nesting property, so that the shape 
of its contour lines is similar to that for t±2 = except 
for the rather rapid decrease for t± > 0.25t. This indi- 
cates that dSC is insensitive to the nesting property (or 
the low-energy property of the spin fluctuation).^' ^'-^^ As 
a result, the superconducting dome emerges. In such a sit- 
uation, the superconducting correlation keeps increasing 
by pressure even when SDW correlation is suppressed, 
as is observed in (TMTTF)2SbF6. The superconducting 
correlation also has a maximum at around t± ^ 0.25t. 
The rather rapid decrease of Xd for t± > 0.25t is due to 
the strongly reducied spin-fluctuations besides the indi- 
vidual interchain hopping. 

Then, we show another case, t±2/t = 109.6 x {t±/t)'^ 
in Fig. 5 (c) , where the Fermi surface keeps good nesting 
condition. The behavior of SDW changes, but the essence 
does not change, namely, the contour lines fall down be- 
low T = t±2- The behavior of dSC is not affected so 
much. 

These contour plots of Xii{T,t±) can be roughly con- 
sidered as P-T phase diagrams by introducing a small 
interplane coupling, V^. In this case when Xi^{T,t±) 
reaches the value V~^, a phase transition occurs, so that 
Fig. 5 represents to the phase diagram. Wo assume that 
the three dimensional transition temperatures are deter- 
mined by the condition 1 — V^x^(T, t±) = with a con- 
stant Vf^. Note that, in order to estimate the real Tc and 
obtain the real P-T phase diagram, we should take into 
account the details of the compounds and the pressure 
effect, such as the rcniornialization of the interplane pro- 
cess to V^, and the pressure effects on U/t and V^. But 
the important statement here is that Xd (7c) essentially 
increases by t± for small which have not been ob- 
tained by a conventional superconducting mechanism for 
higher dimensions, even with a simple assumption of the 
constant V^. Whether Tc kccip increasing after the SDW 
phase is suppressed is determined only by the property 
of Xs, and hence by the energy scale of the nesting de- 
viation T = t±2- Therefore, even with a small diffrence 
of the nesting conditions, the qualitative property of the 
phase diagram can change. When the nesting is easily 
broken by the pressure, the phase diagram is like Fig. 5 
(b). The results of (TMTTF)2SbF6 corresponds to this 
situation. When the nesting property is solid against the 
pressure, the phase diagram becomes like Fig. 5 (c). The 
TMTSF-salts would be in this situation. 

One may notice that the situation of Fig. 5 (b) re- 
sembles that of high-Tc superconductors. Actually, the 



effects of the quantum interference between the Cooper 
and Peierls channels, which leads to the increase of Tc in 
the present case, are also discussed in the context of the 
origin of the pseudo-gap behavior of high-Tc supercon- 
ductors. ^° But for the isotropic square-lattice, the RG 
method yields the monotonically decrease of Tc.^^'^^'^* 
Therefore, in order to understand the Tc behavior of 
high-Tc superconductors, it would be necessary to take 
account of the strong correlation in the vicinity of the 
Mott insulator. ^^"^•^ The present results of Tc is not rel- 
evant to the high-Tc superconductors; it is, instead, a 
characteristic behavior of QID systems. We also note 
that the increase of Tc is also realized when the long- 
range Coulomb interaction is finite and the symmetry of 
the superconductivity is /-wave triplet. 

In conclusion, we have calculated the response func- 
tions both for SDW, Xs, and d-wave singlet supercon- 
ductivity (dSC), Xd, by the A/'-chain RG technique, and 
investigated the effect of dimensionality upon the spin- 
and superconducting-fluctuation. It is shown that, for 
the perfect nesting case, not only Xs but also Xd are en- 
hanced by t±. This increase of Xd is realized only in the 
system where following two conditions are satisfied; i) 
the interference between Cooper and Peierls channels in 
higher energy region, and ii) fc_L-dependence in couplings 
in lower energy region. These conditions are satisfied only 
in QID systems, and not in two (square) or three dimen- 
sional systems. The system gains largest superconduct- 
ing correlations at around t± ~ 0.25, where both ID 
fluctuation and spin fluctuation contribute to supercon- 
ductivity. When the Fermi surface deviates from perfect 
nesting, Xs(T) is cut off at T = t±2, whereas Xd{T) is not 
affected so much. Consequently, Xd increases even after 
the spin fluctuation is suppressed, which leads to the fact 
that Tc keeps increasing away from the SDW phase. This 
is also valid for the /-wave triplet pairings. Our results 
give a possible mechanism of the increasing- Tc behavior 
observed in (TMTTF)2SbF6. 
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